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By using the perturbation expansion up to the fifth order, we study the two-triplet-dimer 
excitation spectra in the Shastry-Sutherland model, where the localized nature of a triplet-dimer, 
the propagation of a triplet-dimer pair by the correlated hopping and the long-range interactions 
between triplet-dimers play an essential role. It is found that the dispersion relations for first- 
neighbor triplet-dimer pair excitations with S = 1 and p-type symmetry qualitatively explain 
the second- lowest branch observed in the neutron inelastic scattering experiment. It is also 
predicted that the second-lowest branch consists of two components, Px- and py-states, with 
slightly different excitation energies. The origin of the singlet mode at 3.7meV observed in the 
Raman scattering experiment is also discussed. 

KEYWORDS: SrCu2(B03)2, Shastry-Sutherland model, two dimensions, bound state, Heisenberg antlferromag 
net, perturbation expansion, dispersion relation 
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There has been a growing interest in low-dimensional 
quantum spin systems, since one can observe a variety 
of properties where classical pictures break down, fn 
the last decade, spin gapped ground states, particularly, 
in two-dimensional systems have received considerable 
attention in connection with the high-Tc superconduc- 
tivity. For instance, a quasi-two-dimensional compound 
CaV409 is known to have a^nin gap originating in the 
plaquette RVB mechanism .BH'Era) 

Two years ago, another new two-dimensional spift,gap 
system SrCu2 (603)2 was found by Kageyama et alB> In 
this compound, magnetic ions Cu^^ {S — 1/2) are ar- 
ranged as shown in Fig. 1. Miyahara and UedaQ^ pointed 
out that SrCu2(B03)2 is an experimental realization of 
a special class of the Heisenberg antifeo-omagnets which 
is called the Shastry-Sutherland modefl^ 



n 



(1) 



where denotes a (next-)nearest-neighbor 

pair of spins. The strength jaLthe interdimer coupling 
is now considered A ^ O.GSBEP The orthogonal dimer 
structure as seen in Fig. 1 leads to the following unique 
properties: the direct product of the singlet-dimer states 
of the nearest-neighbor pairs of spins is the ground state 
exactly for A <~ 0.7 and a triplet-dimer in the singlet 
sea is almost localized.&O 

The nature of the quantum phase transition has 
been qrie of irririnrtant subjects in theoretical investiga- 
tions.&Q'E3Eil'llj-' At first stage of the investigations, it 
was assumed that the exact dimer singlet state is destabi- 
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Fig. 1. Structure of a two-dimensional network formed by Cu'^"'" 
in SrCu2(B03)2. The closed circles represent copper ions. The 
nearest-neighbor bonds are expressed by the bold lines and the 
next nearest-neighbor bonds by the gray lines. The vectors, ex 
and Ey, are the primitive translation vectors in this lattice. Note 
that there exist two types of third-neighbor triplet-dimer pairs, 
which are indicated by 3n and 3n'. For brevity, we call the 
farmer the third-neighbor pair. Within the fifth-order calcu- 
lation, the diagonal interactions between triplet-dimers appear 
for the first-neighbor (In), second-neighbor (2n), third-neighbor 
(3n) and fourth-neighbor (4n) triplet-dimer pairs. 



lized against the Neel ordered state when A is increased. 
A recent study, however, suggested that there exists an 
intermediate phase, where a plaquette RVB state is sta- 
bilized, between the|-£pcact dimer singlet phase and the 
Neel ordered phase.t3 The first-order transition point 
between the exact dimer singlet phase and the plaquette 
RVB phase was estimated as Ac = 0.677 by using the 
series expansion method. Another interesting topics in 
SrCu2 (603)2 is that the magnetization plateaus are-t, 
served at 1/3, 1/4 and 1/8 of the full Cu moment I 
The localized nature of a triplet-dimer plays an impor- 
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tant role in the appearance of the plateaus. HEi 

As mentioned above, a triplet-dimer is hard to prop- 
agate in the singlet sea because of the orthogonal dimer 
structure. In precise, the orthogonal dimer structure pro- 
hibits the propagation o£ a_triplet-dimer up to the fifth- 
order perturbation in A.HEj^ Recently, almost flat band 
of the single-triplet-dimer excitations at S.OmeV was ob- 
served in the neufaon inelastic scattering experiment by 
Kageyama et alS3 

In this neutron inelastic scattering experiment, 
multiple-triplet-dimer excitation spectra were also ob- 
served, which show rather dispersive behavior compared 
with the single-triplet-dimer excitation. In the ESR mea- 
surement, the second spin-gap was observed at 4.7meV, 
which corresponds, to the band bottom of the second- 
lowest spectrum.llS Note that the localized nature of 
a triplet-dimer leads to almost flat bands of multiple- 
triplet-dimer continua and makes the experimental ob- 
servation of various multiple-triplet-dimer bound states 
possible. The observation of a singlet bound state at 
3.7meV in th£_.Raman scattering experiment has been 
also reported.Ej' 

In this letter, we study the two-triplet-dimer excita- 
tions by the fifth-order perturbation expansion. We here 
mention the relation between the present problem and 
the two-magnon bound state in the IsingJEJeisenberg fer- 
romagnets which is now well understood.E3' In the Ising- 
Heisenberg model, when the hopping of a magnon is sup- 
pressed by the Ising-like exchange anisotropy, the forma- 
tion of the two-magnon bound states is more favored and 
the band of those states becomes more narrow. Within 
the first-order perturbation, the mechanism of the for- 
mation of two-triplet-dimcr bound states in the present 
model is the same as that of the the two-magnon bound 
states in the Ising limit of the Ising-Heisenberg model. 

In the Shastry-Suthcrland model, the two-triplet- 
dimer bound state problem becomes to be more interest- 
ing due to the following two reasons arising from higher- 
order perturbation processes. First, the range of the 
diagonal interactions between triplet-dimers is enlarged 
when the order of the perturbation is extended. Sec- 
ond, the correlated hopping, where a triplet dimer can 
hop when the rest spectator triplet-dimer always lies at 
a neighboring site before and after the hopping, is possi- 
ble in the lower-order perturbation processes than the 
six cplpc. jjihich leads to the single-triplet-dimer hop- 
ping.ESEj'ES^ The former separates a branch of two- 
triplct-dimer bound states from the two-triplet-dimer 
continuum. The latter make such the branch a dispersive 
one. It is also important to observe these higher-order 
effects pronouncedly that SrCu2 (603)2 locates near the 
first-order transition point. We expect that the two facts 
mentioned above and the localized nature of a triplet- 
dimer are properly taken into account in our fifth-order 
perturbation calculation. 

We now turn to the calculation of matrix elements in 
the fifth-order effective Hamiltonian, , for the sub- 
space composed of two triplet-dimers. We write the two- 



triplet-dimer states as 
1 

\r,S,S)^ J2 CiS,m)\tm)r\t-,n)r+s]X\')r', (2) 
m= — 1 r' 

where S = 0,1,2. In this definition, \s)^ represents the 
singlet state of the dimer at r, and \tm)r represents the 
triplet state with the total of to. The prime attached 
on the product symbol means the exception of the sites 
r and r + S. The explicit expression of C{S, m) is given 
by 



C = 



-l/\/3 l/\/3 -1/V3 

-1/V2 
i/Ve 2/V6 1/V6 



(3) 



where the row (column) index runs S = 0,1,2 (to — 
— 1,0,1). Note that there is a symmetry relation 
C{S, —to) = (— 1)'^C(5', to). The relative vector between 
the two triplet-dimers, 6, is assumed to satisfy i5j; > or 
Sx — 0, (5y > in the practical calculations. 

The interactions between two triplet-dimers within 
fifth-order approximation are given by 



Vr^siS) = (r, d, S\ |r, S, S) - {Eg + 2A 



(4) 



where Eg — — 3A^d/4 is the ground state energy of the 
system with Nd dimers and 



A,g = l 



,2 A3 



5A5 
"32 



9r (5) 



is the fifth-order series of the spin gap.EfLi' Nonzero el- 
ements of Vj.,s{S) are summarized in Table I. The in- 
teractions for 5 = 2 have been calculated already up to 
the third order and the first-order repulsion between the 
first-neighhaiilriplet-dimers is known to lead to the 1/2- 
plateau.By'EJ On the other hand, we find in Table I 
that the first-order interactions between first-neighbor 
triplet-dimers for 5 = and 1 are attractive. 



Table I. Nonzero elements of the diagonal interactions between 
two triplet-dimers, ^('S')- A set of {r,d) for first-, second-, 
third- or fourth- neighbor triplet-dimer pairs is, respectively, de- 
noted by In, 2n, 3n or An. (See Fig. 1.) 





S 




In 




1 

2 


A(-16 + 8A + 16A2 + ISA^* - A*)/16 
A(-8 + 16A + 14A2 - 9A3 - 16A4)/16 
A(8 + 8A - 2A2 - 9A3 - X'^)/16 


2n 



1 
2 


A=^(-16 + 45A2)/32 
A3{-4-6A + 3A2)/16 
A3{2 4-3A-|-3A2)/8 


3n 



1 
2 


A2(-32 - 48A + 56A2 + 289A3)/32 
A2(-32 - 48A + 8A2 + 129A3)/64 
A2(32 + 48A + 8A2 + 129A3)/64 


4n 



1 
2 


A*(-8 - 17A)/32 
A4(-8 - 17A)/64 
A4(8 + 17A)/64 
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We introduce the Fourier transformation of eq. 



\q,d,S, a) 




J2e"'''-K,6,S), (6) 



where a{= A, B) denotes the sublattice index and Tq, 
denotes a site on the sublattice a. The representa- 
tion of using this basis set becomes to be a direct 
sum of finite-size matrices because of the prohibition of 
the single-triplet-dimer hopping. In the fifth-order ap- 
proximation, the representation of H'^^ has off-diagonal 
matrix elements among first-, second-, fourth- neighbor 
triplet-dimer pair states. For the third-neighbor triplet- 
dimer pairs, has only a diagonal element, which is 
different from Eg + 2Asg. For the other states, H'^^ has 
only the diagonal element of i?g -I- 2 Agg . 

We here mention the zeroth-order wave-functions for 
the first-neighbor triplet-dimer pair excitations, which 
are relevant to the observed multiple-triplet-dimer exci- 
tations in the neutron inelastic scattering experiment as 
described below. The first-order perturbation removes 
the degeneracy between the first-neighbor triplet-dimer 
pair states and the other states. We write the eigenfunc- 
tions of the eigenvalue problem to determine the second- 
order energies as follows: 



V2 



til 
V2 



\q,e^,S,a^) , (7a) 



\q,e^,S,a^) (7b) 



V2 V2 

for ^ ^ x,y, where we have defined Ux = A, Uy = B, 
Oix = B and Uy = A. Note that |s{)5g and 1^5)50 have, 
respectively, the s-type and p^-type symmetries. The 
degeneracy between the s-wave states and the p-wave 
states is removed by the second-order perturbation. The 
fourth-order perturbation removes the remaining two- 
fold degeneracy for each of the s-wave states and the 
p-wave states at the general points in the first Brillouin 
zone. Then the way of mixing between the two s-wave 
functions is determined as follows: 

l^ + )s,<7 = f{<lx,qy) \Sx)s^q + {-'^ffi^y, Ix) \Sy) g^q : 



fiqy,qx) \Sx] 



S:<7 



-^ff{qx,qy)\Sy] 



where 

[f{qx,qy)]-^=l + 

3S qx — cos qy — y^l + (cos qx — cos qyY 



(8b) 



(9) 



The results of the fifth-order perturbation expansion of 
excitation energies, A_B, for all the two-triplet-dimer ex- 
citations are shown in Fig. 2, where we choose A = 0.55. 
In this figure, the results except for the second- and 
fourth-neighbor triplet-dimer pairs are obtained by the 
usual series expansion. For the second- and fourth- 



neighbor triplet-dimer pairs, the mixing among these 
states at high-symmetry points of q occurs in the fifth 
order, which makes the convergence of the series slow 
around such the high-symmetry points. To avoid this, we 
diagonalize directly the effective Hamiltonian for the sub- 
space of the second- and fourth-neighbor triplet-dimer 
pair states, which is an hermitian matrix obtained by 
making further partial diagonalization for . The re- 
sultant eigenvalues are plotted in Fig. 2 as the solid lines. 
Roughly speaking, almost dispersionless bands indicated 
by "4n" originate in the fourth-neighbor triplet-dimer 
pairs. The other solid lines originate in the second- 
neighbor triplet-dimer pairs. 




0.0 
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Fig. 2. Dispersion relations of the two-triplet-dimer states for (a) 
5 = 0, (b) 5 = 1 and (c) S = 2, where we choose A = 0.55. For 
S = 1, dispersion relation of the single-triplet-dimer states is 
also shown. 



We find dispersive behavior for the branches of the 
first-neighbor triplet-dimer pairs for S" = 0, I and 2, 
which can propagate in the singlet sea from the third or- 
der. It should be stressed that the dispersion relation of 
the first-neighbor triplet-dimer pairs with p-type symme- 
try for S = 1 qualitatively explains the behavior of the 
second-lowest branch in the neutron scattering experi- 
ment. (See Fig. 3 in Rcf. 17.) Very recently, a neutron 
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inelastic scattering experiment under the magnetic fields 
was carried out.clP This experiment revealed that the to- 
tal spin of the second-lowest branch is one, which is also 
consistent with our result. Note that the fourth-order 
terms remove the two-fold degeneracy between p^.- and 
Pj,-type first-neighbor two-triplet-dimer pairs at general 
points in the first Brillouin zone. The explicit expres- 
sions of the excitation energies are given by 

AEp^{S = 1) = 2 - ^(1 + cosqxCOsqy) 

A** 

— — (35-1-15 cos cos Qy — 2 cos^ qy sin^ q^ ) 
A^ 

— — [274-1- 175 cos qx cos — 8 cos^ qx 

X {l + (6 + cosga; cos qy) sin^ qy}] , (10a) 

AEp^{S = 1) = {qx ^ qy in rhs of cq. (10a)), (10b) 

for the Px- and pj^-states, respectively. As seen in the 
above equations and Fig. 2(b), the observed second- 
lowest branch is composed of two components with 
slightly different energies on the path from the T point 
to the X point and on the path from the X point to the 
M point, which may be detected by careful analysis of 
the experimental data. 

We turn to the second- neighbor triplet-dimer pairs. 
From Fig. 2, we find dispersive behavior for S* = and 
1, but we find rather flat bands for 5 = 2. Note that 
the effective Hamiltonian for S* = 2 is nothing but the 
two-body part in the subspace of the effective Hamil- 
tonian dominating the magnetization process. Rather 
small energy-gain by the correlated hopping for 5 = 2 
seems to be consistent with the appearance of the mag- 
netization plateaus in SrCu2(B03)2- In Ref. 16, where 
the magnetization plateaus are studied by the third-order 
perturbation, it was assumed that the states with first- 
neighbor triplet-dimer pairs can be truncated. We find in 
Fig. 2(c) that the excitation energies of the first-neighbor 
two-triplet-dimer pairs for 5 = 2 are not the lowest lying 
branch even when the fourth- and fifth-order terms are 
taken into account. This finding gives a support to the 
treatment in Ref. 16. 

As for the third-neighbor triplet-dimer pairs, the bind- 
ing energies are not so small, but propagation does not 
occur within the fifth-order perturbation. The fourth- 
neighbor triplet-dimer pairs can propagate from the fifth 
order but the band width is too narrow to see in this 
figure. 

In Fig. 3, we show the minimum two-triplet-dimer ex- 
citation energies as a function of A. For 5=1, the result 
of the exact diagonalization on a finite-size cluster with 
A^D — 10 is also shown. We find that the series of the 
two-triplet-dimer excitation for 5 = and 5 = 2 tends 
to converge within the present calculation. In particu- 
lar, the fifth-order perturbation series results in 3.7meV 
for the minimum energy of the singlet two-triplet-dimer 
excitation for SrCu2(B03)2 with A = 0.63 and the in- 



tradimer exchange constant of 7.3meV. This value agrees 
with the observed singlet excitation energy in the Ra- 
man scattering experiment .E3' As for 5=1, the result 
in Fig. 3(b) indicates that we should extend the order 
of perturbation to obtain the quantitative description of 
SrCu2 (303)2. Note that the exact diagonalization gives 
4.6meV for the minimum two-triplet-dimer excitation en- 
ergy for SrCu2 (603)2, which agrees with the^ results of 
the neutron inelastic scattering experimentO^ and the 
ESR measurement £3-' Our fifth-order result of this exci- 
tation energy is lower than that of the exact diagonal- 
ization. If higher-order terms beyond the fifth order are 
taken into account, then the binding energy is expected 
to be suppressed due to the single-triplet-dimer hopping. 
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Fig. 3. The minimum two-triplet-dimer excitation energies as a 
function of A for (a) S = 0, (b) 5 = 1 and (c) S = 2. The solid, 
dash-dotted, long-dashed and short-dashed lines, respectively, 
represent the fifth-, fourth-, third- and second-order results. The 
closed circles in (b) are the results of the exact diagonalization 
for the finite-size cluster with A^d = 10. 



In summary, we have studied the two-triplet-dimer ex- 
citations in the Shastry-Sutherland model by using the 
fifth-order perturbation. It has been found that the 
second-lowest branch observed in the neutron scatter- 
ing experiment and the singlet state at 3.7meV observed 
in the Raman scattering experiment originate from the 
first-neighbor two-triplet-dimer excitations. It has been 
also pointed out that the second-lowest branch in the 
neutron scattering experiment is composed of two com- 
ponents, Px- and Pj,-states, whose degeneracy is removed 
by the fourth-order perturbation. We expect that such 
the structure is confirmed by the careful analysis of the 
experimental data. 

The author would like to thank Prof. A. Oguchi for 
useful discussions and critical reading of the manuscript. 
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mori in TITPACK Ver. 2. 
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